Perspectives on Intracluster Enrichment and the Stellar Initial 
Mass Function in Elliptical Galaxies 



ABSTRACT 

The amount of metals in the Intracluster Medium (ICM) in rich galaxy clus- 
ters exceeds that expected based on the observed stellar population by a large fac- 
tor. We quantify this discrepancy - which we term the "cluster elemental abun- 
dance paradox" - and investigate the required properties of the ICM-enriching 
population. The necessary enhancement in metal enrichment may, in principle, 
originate in the observed stellar population if a larger fraction of stars in the 
supernova-progenitor mass range form from an initial mass function (IMF) that 
is either bottom-light or top-heavy, with the latter in some conflict with observed 
ICM abundance ratios. Other alternatives that imply more modest revisions to 
the IMF, mass return and remnant fractions, and primordial fraction, posit an 
increase in the fraction of 3 — 8 M stars that explode as SNIa or assume that 
there are more stars than conventionally thought - although the latter implies 
a high star formation efficiency We discuss the feasibility of these various so- 
lutions and the implications for the diversity of star formation, the process of 
elliptical galaxy formation, and the nature of this "hidden" source of ICM metal 
enrichment in light of recent evidence of an elliptical galaxy IMF that, because 
it is skewed to low masses, deepens the paradox. 



1. Context 

The hot plasma that pervades the volume of the most massive galaxy clusters - the 
intracluster medium (ICM) - provides a wealth of diagnostic data on the process of galaxy 
formation in structures formed from the largest primordial density fluctuations to have en- 
tered the nonlinear regime and undergone gravitational collapse. The history and efficiency 
of star formation, and the effects of interactions among galaxies and between galaxies and 
the environment in the form of infall, outflow, and dynamical stripping, are reflected in the 
thermal and chemical properties of the ICM - both in individual systems and in the evolving 
population of clusters. 

In massive (> 10 14 M Q ) clusters, the ICM dominates the baryon inventory and accounts 
for > 10% of the total matter content (e.g., Lagana et al. 2011). The discovery that cluster 
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gas fractions expected to be representative of the universe as a whole exceed fi&aryonQ precip- 
itated a crisis referred to as the "baryon catastrophe" when combined with the assumption 
of a matter-dominated unive rse and evidence for, and the in flationary prediction of, a flat 
universe - i.e. flmatter = 1 (jFabianlll991t IWhite et al.lll993[ ). Of course this paradox was 
resolved by the discovery of dark energy and the concordance cosmology which reconciles a 
flat universe with a reduced Q, ma ue r- The cosmic bary on matter fraction is now accurately 
determined, ^ baryon / 'Vt matter = 0.17 (IJarosik et al.ll201ll ). with a relat ively modest percent age 
of the baryons collapsed into stars, VL stars /VL matter f» 0.0074 — 0.011 (jGallazzi et al.l 120081 ) . 



However a related paradox persists to this day. While a solar nucleosynthetic yield 
(i.e., star formation ultimately producing the amount of metals necessary to enrich one solar 
mass to solar abundances) is sufficient only to enrich baryons to average abundance of < 0.1 
on a universal scale , cluster baryons are enriched in Fe (and other elements) to ~hal f- 
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ies is evidently more efficient than in the field; however, a large discordance - first em- 
phasized following the groundbreaking accuracy and range of abundance measurements 
made with the ASCA X-ray Observatory - remains between cluster metals and the num- 
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MfpJLr ~ 0.0065 - falling short by a factor of 5 or more compared to what is measured 



(ISakuma et al.ll201ll ; ISato et al.ll2012f ). 



To resolve this "cluster elemental abundance paradox" one generally must conclude that 
there was more star formation in clusters than conventionally estimated, and/or that star 
formation in galaxy clusters has an enhanced efficiency of producing supernova progenitors 
and synthesizing metals. "Missing" stars may take the form of low surface brightness in- 
tracluster light (ICL), infe rred both from o bservations (Zaritsky et al. 12004 ; Lin fc Mohr 



20041 ; I Gonzalez et al.l |2007J) and simulations (IPuchweine al 



2010 



Rudick et al.l 1201 if ). Re- 



cent literature includes a large range in calculated star-to-ICM ratios primarily due to di- 
vergent ICL estimates, but also to different assumed stellar mass-to-light ratios. A stellar 
initial mass function (IMF) that is relatively top heavy increases both the ratio of stars 



f2 represents mean densities in units of the critical density that delineates closed and open universes. 
2 We adopt the solar abundance standard of lAsplund et al.1 (|2009l ) where Fe/H= 3.16 x 10~ 5 by number. 
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formed, and of me t als created, t o present-day light. S uch an IMF may be bimodal in nature 



flElbaz et al.l Il995l : iLarsonl Il998l : Moretti et al. 



2003 ), the seco nd mode perhaps associated 



with a distinct pre-enrichm ent population (IBregman et al.ll2010l ) where Population III hyper- 
novae (ILoewensteinl 120011 ) may play a role. Given the con ventional wisdom that most ICM 
metals originate in elliptical galaxies ( lArnaud et al.lll992[ ). this problem clearly connects to 
the fundamental galaxy formation questions of the IMF in ellipticals and the transport of 
material from these galaxies into intergalactic space. 

In this paper we undertake a fresh and comprehensive, though generic, examination 
of the metal inventory in the ICM of rich galaxy clusters. We focus on addressing a single 
well-defined, though multifaceted, question: what characteristics of the stellar population are 
necessary to produce the observed level of ICM enrichment? In doing so we address issues 
related to the IMF, star and galaxy formation efficiency, galactic winds, the astrophysics of 
supernova progenitors and explosions, and the apportionment of products of different super- 
nova types into stars and ICM. Section 2 quantitatively summarizes the cluster elemental 
abundance paradox. 

In Section 3, where a standard IMF is assumed, the level of enrichment and abundance 
pattern are related to phenomenological parameters that encapsulate the star formation 
efficiency, and the demographics of supernovae and the success of stars in locking up the 
products of their explosion. Substantial departures from standard values are required to 
match observations. In Section 4, we cast a wider net by considering the effects on ICM 
enrichment of a wide range of IMFs in the context of a self-consistent galaxy chemical 
evolution treatment that accounts for the relevant astrophysical quantities. When juxtaposed 
with recent evidence for an IMF in elliptical galaxies that produces fewer metals than a local 
IMF, this analysis reinforces and clarifies the conflict between the ICM metallicity and the 
characteristics of the stellar population generally assumed responsible for ICM enrichment. 
Results and their implications are discussed, and conclusions summarized, in Section 5. 



2. The Cluster Elemental Abundance Paradox Quantified 



2.1. Basics 



Consider the total baryon mass in a cluster of galaxies within some sufficiently large 
radius that it may be considered a closed box in the chemical evolution sense - that is, 
all products resulting from the transformation of some of this gas into stars (including the 
stars themselves) are contained within this radius. There is direct evidence that this is a 
good approximation for sufficiently massive clusters if invoked at radii that are a significant 
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fraction of the virial radius, based on the consi stency between the total cluster baryon 
fraction and the universal value mentioned above (ILandry et al.ll2012[ ). perhaps with a 10% 
"depletion" correction at r^nn (the radius within which the ave r age m a ss density is 500 time s 
the critica l dens it y): see iGonzalez et al.l toQl\ \: IPratt et all J2009h : iGiodini et all d2009h : 
Ade et al. ( 2012 ); Planelles et al. ( 20121 ) . Presumably this is a result of the extreme depth 
of their gravitational potential wells. However it is possible that this also applies to galaxy 
groups, and perhaps even giant elliptical galaxies if one could inventory the gas all the way 
out to the virial radius (and perhaps beyond, if entropy injection has dispersed the gas 
distribution). 

We define the overall efficiency of converting gas into stars, e s f, such that the total mass 
in stars formed (regardless of where) is 
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*,form 



sf -^-baryon: 



X 



where M^ aryon is the total baryon mass being considered. At the present time the total mass 
in stars, whether contained in individual cluster galaxies (including the brightest cluster 
galaxy - BCG) or associated with intracluster light (ICL), is 



M^baryon &sf (1 



(2) 



and the mass in gas is 



M, 



gas 



-^■baryon M-* 
■^baryon [1 



r*)l 



(3) 
(4) 



where r* is the stellar "mass return fraction" - the fraction of the mass previously formed into 
stars recycled back into gas. For massive clusters one may neglect the distinction between 
the total mass in gas and the mass in the ICM, and henceforth we equate M gas with Mjcm- 
The star formation efficiency, in terms of the observable Micm/M* is 



£ sf 



:i-r* 



1 + 



M 



ICM 



-1 



(5) 



We consider the enrichment of cluster baryons in chemical elements released in supernova 
explosions, i.e. those of atomic number A > 8. It is straightforward to extend the analysis 
to elements statically synthesized in intermediate mass stars; however, the elements (C, N) 
to which this applies are not well-constrained by current X-ray observations. 

Both the overall level of baryon enrichment (that is, the metallicity) and the abundance 
pattern are determined by the total number of supernovae and their nucleosynthetic yields. 
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We separate the enrichment contributions of the two main classes of supernovae - Type 
la supernovae (SNIa) that result from the explosion of a white dwarf, and core collapse 
supernovae (SNcc). Their total numbers may be expressed as 

N cc = r] cc M*j orm (6) 

and 

N Ia = V Ia M* Jorm , (7) 

where rf c and r/ Ia are, respectively, the specific numbers of SNcc and SNIa explosions per 
star formed. It is useful to define the total supernova number, supernova ratio, and SNIa 
fraction as follows: 

N SN = N cc + N Ia (g) 

= V SN (l-O^M*, (9) 



where 



and 



V SN = n cc + V Ia ; (10) 

la 

R SN = ^-; (11) 

la pSN 

f Ia = V_ = _^ (12) 

J ~ V SN 1 + R SN ' { ' 
Note that the number of SN per unit mass in the ICM is 

^ = ^(1-^7^, (13) 

MlCM J ICM 

where /* and ficu are the present-day mass fractions of stars and gas: /*/ ficu = M*/M IC m- 



2.2. Stars and Supernovae 

A combination of theoretical and empirical considerations enter into determination of 
the stellar and supernovae parameters. Star formation is not sufficiently well-understood to 
allow an a priori estimate of e s f, and M* is estimated from observations of cluster starlight 
in galaxies and intracluster space. To convert to mass, stellar population synthesis can be 
employed, involving assumptions about the IMF - among many other factors. The mass- 
to-light ratio in individual galaxies can also be inferred from dynamical modeling of stellar 
velocity dispersion distributions. The mass return fraction, r*, may be calculated from the 
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IMF, star formation history (SFH), and relation of stellar remnant mass to progenitor mass 
for individual stars estimated from stellar evolution theory and observations in the Galaxy. 
The SNcc efficiency, 7] cc , depends on the IMF and range of masses that result in core collapse 
explosions. Although one may model the time-dependence of the SNIa rate from assumptions 
about the binary star progenitor configuration and the distributions of binary mass ratios 
and separations, the normalization is difficult to estimate a priori. Empirical estimates of r) Ia 
(and 7] cc ) are emerging from supernova surveys that al so constrain the distribut i on of delay 



times from the observed evolution of the SNIa rate (IMaoz fc Mannuccil 120121 ; ISand et al. 
20121 ). Estimation of these fundamental quantities - M*, r*, 77 , rf c - in principle requires a 



convolution and a time-integration over an evolving galaxy population with disparate SFHs 
and, conceivably, IMFs. A first order approximation treats the ensemble stellar population as 
a single simple population of stars formed at some (high) redshift with a common IMF - an 
approximation most suitable for rich clusters where the total galaxy mass is most dominated 
by elliptical galaxies. 

We can insert some reasonable values to get a sense for the expected level of super- 
nova enrichment and relative contribution from the two classes of supernovae. For a "diet 
Salpeter" IMF that represents a simple alteration - proposed a s a means of reconcil iation 
with the observed relative frequency of ^subsolar mass stars flBell fc de Jong l200lh - of 
the c lassic single-slope Salpeter function (ISalpeterl I1955T). n Ia ^ 
20121 ). r* ~ 0.35 for an old stellar population fjFardal et al.l 2007: 
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0.01. For massive clusters 
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with substantial systematic uncertainty and study-to-study variation 


see ] 


^eauthaud et al. 



(120121 )). and evidence of an increase in magnitude and scatter with decreasing cluster 
mass. The resulting total number of supernova explosions per solar mass of ICM is ~ 
1.5 x 10- 3 (10/,// /CM ) M- 1 . 



2.3. Application to the ICM 

The equation for the mass of the ith element in the ICM, Mj, in terms of the number 
of SNIa and SNcc that enrich the ICM, N Ia and N cc , and the yields per SNIa and SNcc, y( a 



-7- 



and (yf), is 

^ = N cc (y?} + N Ia y{ a = N SN (I + R? N Y\(y?) + R SN y! a ), (14) 

where, as defined above, R SN = N Ia /N cc and N SN = N Ia + N cc . The IMF(0)-averaged 
SNcc yield is 

f m " p dmd>(m)yf c (m) 
J m amcpym) 

where m cc and the lower and upper limits for the masses of SNcc progenitors, and a 

single universal set o f SNIa yields is assumed. Despite a plethora of I MF parameteriz ations 
( jKroupa et al.ll2012l ). there is general consensus that a Salpeter slope (jSalpeterl 119551 ) . ~ 
m -2 ' 35 applies at the high mass end relevant for SNcc - at least for star formation under 
"normal" conditions. 

The resulting mass fraction in the ICM (mass Mjcm) of the ith element, fi, is 

= = + ^T^n + R SN vh, (i6) 

M ICM M ICM 

and the mass fraction relative to the solar mass fraction, 

fi _ N SN (1 + RSNrl f{p_ + R SNVt\ (17) 



/iS M IC M \ /is /id 

The relationship between mass fraction, /$, and abundance, Zi, of the ith element (the number 
of atoms of element i relative to that of H, i.e. the entries in standard abundance tables) is 
z i — (/i/^)(A/^-i?) _1 ) where X and Ah are the hydrogen mass fraction and atomic weight 
{Ah = 1.008 AMU) and Ai the atomic weight of the ith element. Relative to solar, the 
abundance is 

Z A = ~ (18) 

z i& X /j fi Q 

where fi & = z^X^Ai/ Ah), and the approximation X = X & is invoked - an approximation 
that is valid as long as the total mass fraction of metals (< 2% for solar abundances) is 
small and th e He abundance is fix ed. X for various solar standard abundance sets is given 
in Table 4 of ksplund etaJ fcpQ9h . 

Finally, the abundance relative to solar is expressed as 

Zi = — = ^-(i + R SN r\(vf) + R SN y{ a ') (19) 

where N SN / M IC m is given by equation (13); and, = (y^ c ) / f iQ and y\ a ' = y\ a j 'f iQ . For 
Mjcm measured in M Q these are the yields of the ith element relative to the mass of that 
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element contained in one M of solar abundance material. As noted above for the total 
baryons, R SN and N fully determine the level and pattern of ICM enrichment - modulo 
se ts of SNIa and SNcc yields - and can be compared to ICM abundances. The new approach 
of iBulbul et al.l ( 120121 ) directly fits X-ray spectra to a model parameterized by R SN and N SN 
via the abundance predicted by equation (19). 

Focusing on Fe, the element with the most widely determined and most accurate global 
ICM abundance measurement, equation (19) predicts Z Fe j C M = 0.255(10/*///cm); f° r 
the values of R SN and N SN derived at the end of the p r evious subsection and adopting 
{yl a , {yf))= (0.743 M , 0.0825 M ) from iKobavashi et all (|2006h - about half the typical 
observed value for f*/ ficM = 0.1. However this assumes that all of the metals produced by 
supernovae reside in the ICM. The values R SN and N relevant here correspond to those 
supernova explosions that enrich the ICM (or, for some particular X-ray measurement, those 
in a particular spectral extraction region of a particular cluster). Not all of the products 
resulting from supernova nucleosynthesis are available to enrich the ICM. 



2.4. Metals Locked Up in Stars 

One approach to evaluating galaxy cluster enrichment is to estimate the total inventory 
of metals in stars and in the ICM in the context of the total required number of supernova 
explosions. However, our focus here will be on the ICM which offers more accurate abundance 
determinations over a wider range of elements via X-ray spectroscopy. This specifically 
requires a correction accounting for how supernova products are apportioned among gas and 
stars. 

The galactic mass in rich clusters is dominated by early-type systems that form their 
stars rapidly. This results in the well-established enhancement in [a/Fe], the abundance 
ratio of a-elements to Fe (expressed as the logarithm with respect to solar) - i.e., SNcc 
products are pref erentially locked up in stars . In their investigation of the giant elliptical 



galaxy NGC 4472, lLoewenstein fe David (120101 ) found that a ratio of SNIa to total supernovae 



of N^/N™ ~ 0.11 (N?/N? N ~ 0.89) and number of supernova per mass in (present-day) 
stars of N? N /M* ~ 0.0083 resulted in [a/Fe]* ~ 0.25 and Zp e ,* ~ 1 (as observed in this 
particular galaxy, but typical of the class; see also Lin et al. 2003, Gallazzi et al. 2008). This 
enables us to estimate the lock-up corrections, rf^ and 7^ a , needed to convert the specific 
supernova numbers to those available for enrichment of the ICM as follows: 

N Ia 

Vl a = ^r(l - r.) « 6.0 x 10~ A Z Fe ^ (20) 
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and 



4.8 x lQ- 3 Z Fe 



The corresponding values available to enrich the ICM, rjj^ M and rff CM , are 



(21) 



(22) 



and 



Vicm 



v cc - V? 



3.2 x 1(T 3 M 



r-l 



(23) 

for the default parameters considered above. That is, the metal production from ~ 60% of 
SNcc and ~ 30% of SNIa must be locked up in stars to enrich them to solar Fe abundances and 
[a/Fe]* ~ 0.25 - with these factors deducted from the total to obtain the effective enrichment 
of the ICM. A relative overabundance of SNIa contributing to ICM enrichment is expected 
based on the inference that the accelerated formation of stars in clusters preferentially locks up 
the products of SNcc . This asymmetry is observed in the abundance patterns in cluster cores 
(Ide Plaa et al.l 120071 ; lLovisari et al.l l201ll ) , although whether this extends globally remains 
an open question - e.g., the smothering of galactic winds in central dominant galaxies may 
skew the pattern, relative to the ICM as a whole, via concentrated direct injection of SNIa. 

Based on these estimates, for the supernovae remaining available to enrich the ICM, 
RSN „ 0.44 _ Q .30) and N SN /M ICM ~ 0.71(10 f*/f WM ) x 10~ 3 M© -1 . This enriches 
the ICM in Fe to the level Z Fe j C M = 0.155(10/*///c , m)) thus quantifying the paradox that 
baryons in clusters of galaxies are enriched beyond what is expected based on the stars we 
see in galaxies today - unless either the star formation efficiency exceeds that in the field by 
a factor of ~ 3, or supernovae are more efficiently produced per unit star formation. 



3. A More Comprehensive Examination (I) 

In this, and subsequent, sections we investigate stellar and ICM abundance predictions 
for a range of elements - focusing on a subset selected on the basis of a combination of 
accessibility and diagnostic power: O, Mg, Si, Fe, and Ni. First, we consider a wide range 
of published yield sets and apportionment of metals into stars and ICM in an effort to place 
robust constraints on the required efficiency of star formation. In order to be as general 
and assumption-free as possible we introduce two parameters that gauge the efficiency with 
which stars may lock up supernova products and do so asymmetrically, i.e. preferentially 
for SNcc relative to SNIa. We consider the specific effects of varying the IMF, with its 
coupled impact on r*, r] cc , and r] Ia , in the following section, fixing these parameters at the 
values described above (r* ~ 0.35, r] Ia ~ 0.002, rf c ~ 0.008) for immediate purposes. Note 
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that with these parameters, supernovae enrich cluster baryons in Fe to a mass-averaged 
abundance of Z bar ^ Fe = 1.66e s f. 

We generalize our treatment of quantifying the fraction of metals synthesized by super- 
novae that are inaccessible for ICM enrichment due to lock-up in stars by defining a SNcc 
lock-up fraction, /3 CC , and supernova asymmetry parameter, a SN , such that 



P cc =\, (24) 
r] cc 



..SN I n 



and 

a ™ = {w*) ' (25) 

where Rf = r)l a /r)l c . By definition, (3 CC < 1, and a SN > 1 is expected under conditions 
where rapid conversion of gas into stars results in preferential incorporation into stars of 
SNcc products with respect to those from SNIa products that are released over a relatively 
extended time interval. We provisionally adopt the values that correspond to the estimates 
of the previous section as standard for the remainder of this section: /3 CC = 0.6 and a SN = 
2. ICM abundances are calculated, for a given star formation efficiency e s f, from ICM- 
specific versions of equations (5), (10), (11), (13), (19), with the supernovae per star formed 
effectively reduced to 

Vicm = 7** (l - ^) , (26) 

and 

Vwm = V cc (l - n- (27) 

The results, assuming Kobayashi et al. ( 2006 ) supernova nucleosynthetic yields J§ are dis- 



played in Figure 1. Figures la and lb confirm that the default (IMF and) lock-up parameters 
predict a stellar population enriched to solar Fe abundances, and ICM Fe abundances lower 
than observed for e s f ~ 0.15 (f*/ fwM — 0.11). One may recover the observed level of ICM 
Fe enrichment for e s f ~ 0.3 (/*/ ficM = 0.24), but only for extreme values of the lock-up 
parameters that imply that most of the Fe produced by stars ends up in the ICM, and that 
stellar Fe abundances are well below solar - in contradiction to observations. 

Figure lc plots Z Fe j C M and Z Fe ^ versus e s f for selected pairs (/3 CC , a SN ) - including 
the default. Also shown is the limiting case /3 CC — > where 100% of metals produced by stars 
reside in the ICM. In this case 

1.66e s f 

J Fe,ICM = : — — — , (28) 



r*) 



3 Although these are averages for an IMF with a Salpeter slope at the high mass end, we adopt these in 
general - the differences for the IMFs we consider are generally small. 
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Fig. 1 . — Left panel (a) : ICM (solid lines) and stellar (broken line) Fe abundance versus lock-up fraction 
(3 CC for (the adopted standard) a SN — 2, and £ s t — 0.15 or 0.30. The horizontal dotted line shows the 
typical observed value abundance, Zp e ,iCM — 0.5, the vertical dotted line the adopted standard (3 CC = 0.6. 
Middle panel (b): Same as (a) versus a SN for (the adopted standard) (3 CC = 0.6. The vertical dotted line 
shows the adopted standard a SN = 2. Right panel (c): ICM (curves) and stellar (horizontal lines) Fe 
abundance versus e s f for pairs (a , /3 cc )=(10,0.1)(short-dashed), (3,0.2)(dotted), (2,0.6) (standard: solid), 
(2.0,0.3)(long-dashed), and (1,0. 8) (dot-short-dashed). The dot-long-dashed line shows the limiting case 
/3 CC -> 0. 



from which on can see that e s j > 0.25 represents an absolute lower limit to the star formation 
efficiency required to enrich the ICM to Z Fe j C M > 0.5. This figure provides an alternative 
demonstration that, for e s f = 0.15, Z Fe j C M < 0.3 even for extreme models where such a 
large fraction of supernova-produced metals is released into the ICM that insufficient metals 
remain available to enrich the stars to the observed level. Both relatively large star formation, 
and small lock-up efficiencies^ are required to simultaneously enrich the stars and ICM to 
the observed level (see, also, Sivanandam et al. 2009). 

The increasing divergence of stellar and ICM abundance ratios (that are independent of 
e s f) with increasing a SN is shown for f3 cc = 0.6 in Figure 2a. One can see how, in this case, 
the enhanced [a/Fe] measured in the old stellar populations that dominate cluster galaxies 
implies a large asymmetry parameter and, as a result, subsolar ratios of a-elements with 
respect to Fe in the ICM. Since [a/Fe]icM « 1 is not observed, large values of a SN may be 
ruled out. This divergence narrows with decreasing (3 CC and ~solar ICM abundance ratios 
emerge for f3 cc ~ 0.3, i.e a lower lock-up fraction (Figure 2b). 

The effects on cluster enrichment of adopting different S NIa yield sets is shown in Figures 



3-5. The solid lines correspond to previous plots that utilize iKobayashi et al.l (120061) yields. 



the dotted lines use the same SNcc yields, but alternative SNIa yield sets from lNomoto et al. 



4 or, more precisely for Fe, f3 cc /a SN << 1 
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Fig. 2. — Left panel (a): ICM (solid lines) and stellar (broken lines) abundance ratios with respect to Fe 
for O (red), Mg (blue), Si (green), and Ni (yellow) versus a for the standard f3 cc = 0.6. Right panel (b): 
same as (a) for (3 CC = 0.3. 



(1997 


); 


Maeda et al. 


(2010 


) - see 



Si (O and Mg are always dominated by SNcc enrichment and insensitive to this choice). 
Figure 3 shows the total baryon enrichment, which is independent of /3 CC and a SN , and 
demonstrates that star formation efficiencies e s f ~ 0.25 — 0.5 are required to enrich cluster 
baryons to a relatively modest Fe abundance of half-solar - e s f ~ 0.4 — 0.7 to attain 0.75 
solar. Figures 4 and 5 that, respectively, show the ICM abundances for the standard lock-up 
parameters, and the maximum ICM abundances corresponding to zero stellar metallicity, 
demonstrate that the conclusions about the required efficiency o f star formation are not a 



result of a particular choice of yields sets. Using SNcc yields from lWoosley fe Weaver! (119951 ) 
(their standard explosion energy, solar abundance progenitor model) does not alter these 
conclusions. 

From the results in this section, we confirm that ICM Fe abundances cannot be produced 
if e s f ~ 0.15 (/*/ ficM ~ 0.11) and quantified the shortfall as a function of how efficiently 
metals in general, and SNIa products in particular, are locked up in stars. For e s f ~ 0.3 
(/*/ fiCM ~ 0.24), they can - but only for small lock-up fractions such that ~ 85% of 
SNIa, and ~ 70% of SNcc, metal production is embedded in the ICM and (by implication) 
ZFe,* ~ 0.5. The assumption of a standard IMF is adopted throughout this section, an 
assumption we relax in the following section. 
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Fig. 3. — Left panel (a): Average Fe abundance of all cluster baryons (stars and gas) as a function of star 
formation efficiency, e s f, assuming a single SNcc yield set but a wide range of SNIa yield sets (see text for 
details). The horizontal dotted line shows the very conservative value of Zp e ^ ar — 0.5. Middle panel (b): 
Same as (a) for Ni. Right panel (c): Same as (a) for Si. 




Fig. 4. — Same as Figure 3 for ICM abundance, assuming the standard /3 CC = 0.6, a SN = 2. 

4. A More Comprehensive Examination (II): Effects of Changing the IMF 

The initial mass function (IMF) of stars in cluster galaxies and intracluster space is 
intimately connected to estimates of ICM enrichment, impacting the mass in stars calculated 
from the total light, the mass return fraction (or, equivalently, the ratio of current stellar 
mass to mass converted into stars), and the numbers of SNcc and SNIa explosions expected 
per mass formed into stars. The general characteristics of the IMF in va rious Milky Way sub - 
populations are now well-determined and generally mutually consistent ( Bastian et al. 201oh. 



Altho ugh several functional forms are commonly used for the "canonical" IMF (IKroupa et al. 



20121 ). these must share the properties of a Salpeter-like slope at high mass with a break to a 
flatter slope below ~ 0.5 — 1 M . Many of the best-studied environments are consistent with 
the hypothesis that this form is "universal" in space and time, but variations in extreme 
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Fig. 5. — Same as Figure 3 for the maximum ICM abundance, i.e. (i cc = 0. 

star forming environments that predominate in the early universe (when most stars in the 
elliptical galaxies that dominate the stellar content in clusters form) may explain a number of 
anomalies, such as an apparent inconsistency between the observed evolution of the global 
star formation rate and stellar mass densities (Narayanan & Dave 2012a, and references 
therein). 

Several recent observational investigations of elliptical galaxies find direct evidence in 
the mass-to-light ratio for either an excess of stellar remnants as realized in a "top-heavy" 
IMF; or, of low mass stars as realized for a "bottom-heavy" IMF (e.g., Cappellari et al. 
2012a; see Section 5 below). By exploiting the level and pattern of ICM abundances we 
constrain the properties of the enriching stellar population. We may then exclude particular 
elliptical galaxy IMFs under the parsimonious assumption that these optically studied stars 
originate from the same parent IMF as those that enrich the ICM - or, alternatively, call 
this assumption into question. 



4.1. Models and Parameters 



Simply put, the level of metal enrichment of the stellar and ICM baryonic sub-components 
in clusters is a reflection of their respective total masses, the total numbers of SNla and SNcc 
that enrich each constituent, and the nucleosynthetic yields of each of these supernova explo- 
sion types. In previous sections, these are expressed in terms of the mass return fraction, r*, 
and formation efficiency, e s f of the stars, the total and relative numbers of SNla and SNcc 
per star formed (R SN , f Ia ), and phenomenological supernova lock-up and asymmetry pa- 
rameters ((5 CC and a SN ) characterizing the ultimate destination (ICM or stars) of supernova 
products. In Appendix A these are further deconstructed into more fundamental astro- 
physical functions and parameters directly connected to stellar and galaxy evolution, thus 
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providing a self-consistent astrophysical framework for understanding the effect of varying 
the IMF on cluster enrichment. Ultimately these are reduced to the following: (1) The func- 
tional form of the IMF (see below); (2) the present-day main sequence turnoff mass (0.9 M ); 
(3) the remnant-progenitor mass relationships for (< 8 M ) intermediate- (equation A4) and 
(> 8 M ) high-mass (equation A29) stars derived, for the former, from well-established stan- 
dard white dwarf masses and, for the latter, from a model for the evolution of the stellar 



population and the delayed-explosion compact remnant prescription in iFryer et al.l ( 120121 ); 
(4) the ratio of mass ejected from galaxies into the ICM during star formation to the mass 
of stars formed, Sow'-, (5) the galaxy formation efficiency, e ga i]o and, (5) various supernovae 
switches and parameters that we now describe. For SNcc we assume progenitor masses from 
m cc = 8 M Q to m up , where m up may differ from the IMF upper mass limit mu (but is the 
same by default, and assumed so in ca lculating the hig h mass return fraction). Since we 



adopt the IMF-averaged SNcc yields of Kobayashi et al.l (120061 ) as a function of progenitor 



metallicity Z CC) Z cc must be specified as well. For SNIa we consider the yield sets described 
in Section 3, assume progenitors in the 3 — 8 M range, and must specify the efficiency e Ia 
defined as the fraction of 3 — 8 M that result in SNIa. In addition, the "prompt" fraction 
of SNIa that explode during the star formation epoch and so may be incorporated into stars 
or ICM, fp 0- , (while a fraction 1 — strictly enrich the ICM) must be specified. 



Our default IMF is the lKroupa et al.l (120121 ) segmented power-law with slopes and mass 



scales that can explain local star formation (Appendix A); other defaults are Z cc = 1, W7 
yields, e gal = 0.25 and 5 GW = 0.5 (e sf = 0.17, f'Jf IC M ~ 0.11), e Ia = 0.076 (Section A.2), 
and fp 0, = 0.5. Under these conditions, r* = 0.41 (r* m = 0.25, = 0.16 for intermediate- 
and high- mass stars, as delineated above) while the fraction in stellar remnants is 0.17 
(0.11/0.06 from intermediate/high-mass stars), if c = 0.011, i] Ia = 0.0022, a SN = 2 and 
(3 CC = 0.39 - similar to the default parameters in Section 3. The resulting abundances are 
shown in Table 1 where, once again as expected, we find reasonable stellar abundances but 
ICM abundances too low by a factor of ~ 2. For comparison we also display the abundances 
for a Salpeter IMF, which fails to provide sufficient metals for all components (including 
stars) for this default set of parameters!^ 

Our approach to examining the effects of varying the IMF on ICM enrichment, that 
attempts to make comparisons at fixed values of the observables to the extent possible, is 
as follows. As detailed in Appendix A, we consider departures from the "canonical" IMF 



Defined as the fraction of baryons initially in galaxies, this is related to the star formation efficiency 
defined in equation (2) by the expression e s f — £gai(^ + Saw)' 1 - 

6 It should be noted that most observational estimates of f*/ficM are IMF-dependent, and would be 
larger for a Salpeter IMF given a fixed amount of optical light. 
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Table 1. ICM Abundances from Default Parameters 








Mg 


Si 


Fc 


Ni 


Baryons 
Stars 
ICM 


0.44 (0.27) 
1.74 (1.05) 
0.30 (0.15) 


0.30 (0.18) 
1.19 (0.72) 
0.20 (0.10) 


0.38 (0.23) 
1.35 (0.82) 
0.27 (0.15) 


0.32 (0.20) 
0.87 (0.54) 
0.26 (0.15) 


0.84 (0.53) 
1.90 (1.21) 
0.72 (0.43) 



Note. — All abundances relative to jAsplund et al.ll2009(l solar standard. The 
values in parentheses are for a Salpeter IMF with the default values of e ga i , &gw > 
e Ia , and f~ a , and the same range of masses (0.07 — 150 Mq). 



Table 2. Model IMF Parameters 











Chi 


sl-0 


0.07 




1.8 -> 2.3 


= "io 


sl-1 


0.07 


0.5 


0.3 -> 2.4 


2.3 


sl-2 


0.07 


1.0 


0.3 -> 2.3 


2.3 


sl-3 


0.07 


0.5 


1.3 


1.8 -> 2.5 


sl-4 


0.07 


1.0 


1.3 


1.8 -> 2.5 


sl-5 


0.07 


8.0 


1.3 


1.5 -> 3.3 


m-1 


0.01 -> 1.0 


1.0 


1.3 


2.3 


m-2 


0.01 -> 0.5 


0.5 


1.3 


2.3 


m-3 


0.07 


0.5 -> 5.45 


1.3 


2.3 



Note. — mhi = 150 Mq for all models displayed here. 
mi, r - equivalent to ml = 0.5 Mq for models sl-1, sl-3, and 
m-2; and 7713 (default: 1 Mq) for models sl-2, sl~4, sl-5, m-1, 
and m-3 - is defined as the mass where the IMF slope transi- 
tions to its high-mass value. Models sl-4 an( i s ^~5 include an 
additional break from a\ = 1.3 to a.2 = 2.3 at m2- eti = ct; 
and 03 = cthi, while ct2 = O( for models sl-2, m-1, and 
m-3; 02 = 03 for models sl-1, sl-3, m-2; and, is set at the 
default C(2 = 2.3 for models sl-1, sl-4, sl-5, and m-2. Param- 
eter ranges correspond to those with physical solutions; i.e., 
fp a and 5qw > & n d e gal < 1 ~ f° r the sl-4 model with 
higher SNIa efficiency (star-to-gas ratio) the range shifts to 
1.5 — ► 2.5 (2 2.5); see below and Tables 3 and 4. 
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(Table 2) described according to sequences with either a single slope below the break mass 
m 3 , or an additional distinct slope below 0.5 M Q . Either the lower mass limit (mi ), or the 
slopes (oti or a 3 ) at low or high mass ends may be varied. Additionally, m 3 may vary in 
sequences of the first type. To isolate the effect of the IMF on ICM abundances, we impose 
invariance on the stellar Fe and O abundances (Table 1). This determines the necessary 
adjustments in the parameters 5gw and f^ a (for some e Ia ). For fixed yield sets, invariance 
in Zp e ,* an d -Zo,* assures invariance in all stellar abundances. Finally, we consider these 
variations at fixed present-day baryon inventory (11% stars, 89% ICM), which is equivalent 
to adjusting e ga i so as to maintain constant £ s /(l — r*) (equations 2, A26) - thus enabling 
us to investigate what adjustments in IMF (if any) may explain observed ICM abundances 
assuming this nominal star-to-gas ratio. The imposition of these constraints rule out those 
IMFs that imply unphysical values of 5gw(< 0), fp a (> 1), or e ga i(> 1) - i.e., some IMFs are 
incompatible with observed stellar abundances and a ~9:1 ratio of ICM to stars (thus the 
limited range of the variable IMF parameters in Table 2). 



4.2. Impact of Varying the IMF 

We remind the reader that our standard IMF has slope a = 1.3 below 0.5 M , and 
2.3 above. Figures 6a-d show the impact on ICM enrichment of varying one (and only one) 
of the slopes and (in some cases) adjusting the break mass, by plotting the Fe abundance 
and Mg/Fe, Si/Fe, and Ni/Fe ratios versus the deviation in the non-fixed slope from these 
standard valuesH This covers many of the IMFs considered in the literature as possibly 
resolving various conflicts between expectations and observations of stellar populations in 
elliptical and/or starburst galaxies. For models sl-1 (si- 2) a below 0.5 M (1 M ) is varied 
with a above these single break masses fixed at 2.3 - i.e. positive (negative) Aa corresponds 
to bottom-heavy (-light) IMFs. For models si- 3, sl-4, and si- 5, a is varied above break 
masses of 0.5 M Q , 1 M , and 8 M , respectively, with a = 2.3 between 0.5 M and the 
break mass in the latter two. For these models, positive (negative) Aa corresponds to top- 
light (-heavy) IMFs. In addition, we plot the results for a single-slope IMF (sl-0). We can 
see that Zpe,iCM ~ 0.5 is predicted for either a (1) bottom-light IMF with a(< 1M ) ~ 1, 
(2) top-heavy IMF with a(> O.5M ) ~ 2 or a(> 1M ) ~ 1.8, (3) single-slope IMF (i.e. both 
bottom- and top-heavy) with a ~ 1.8. As expected, Zpejcu ~ 0.1 for the pure Salpeter 
IMF (single slope and bottom-light cases with slope 2.35). 



7 Mg is almost exclusively synthesized in SNcc, Si primarily (but not exclusively) in SNcc, Fe in both 
SNcc and SNIa, and Ni primarily in SNIa. 
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Varying the lower mass cutoff provides an alternative means of producing either bottom- 
heavy (mio < 0.07 M Q ) or bottom-light (m/ p > 0.07 Mm) IMFs , while increasing the break 
mass that delineates a = 1.3 from a = 2.3 ([Narayanan &: Pavel l2012al ) is also bottom light 
in the sense that the fraction of low-mass stars relative to those at intermediate and high 
mass is suppressed. As shown in Figure 7a, these alternatives can explain Zp e ,icM = 0.5 for 
IMFs with the standard high mass slope a = 2.3 if the IMF lower mass cutoff is shifted from 
the standard mi = 0.07 M Q to ~ 0.2 M Q for an IMF with a single break at 1 M Q , what 
is (essentially) a Salpeter IMF with mi ~ 0.5 M , or for an IMF with the break between 
a = 1.3 and a = 2.3 shifted from 0.5 to ~ 1 M . The last would seem to represent a 
particularly modest departure from the standard IMF. 

Bottom-light and top-heavy scenarios may be directly distinguished in ICM spectra via 
abundance ratios, as demonstrated in Figures 6b-d, 7b- d, and 8a-c; and Table 3. These 
essentially define two branches in the Zpejcu — (ot/Fe)icM plane (Figure 8), with ratios 
in the top-heavy branch connecting to the stellar ratio as the IMF flattens (and — > 1), 
and abundances for the bottom-light IMFs (assured to have f^ a = 0.5) universally rising 
in lockstep with the elimination of low-mass stars. Table 3, confined to those models that 
predict Z Fe j C M = 0.5, illustrates how the a/ 'FepejcM ratios might be exploited to distin- 
guish among bottom-light and top-heavy IMF explanations for ICM enrichment - for the 
former, ratios deviate more strongly fro m those in stars (smaller a/ FeFejcM, larger Ni/Fe) 
and provide a better match to the data (jSimionescu et all 120091 ). 



4.3. Implications of Models with Nonstandard IMFs 

Table 4 displays the essential characteristics of models that produce Z Fe j C M = 0.5, and 
are constrained to match the standard values of stellar metallicity and /*/ '//cmE Several 
general properties, as well as others that distinguish top-heavy from bottom-light solutions 
emerge. Relative to the model with canonical IMF, all have a relative deficiency of unevolved 
low-mass stars, and hence ~ 30% higher mass return (r*) and remnant (/ rem ) fractions - 
implying ~ 30% upward adjustments in the integrated mass of stars formed based on the 
present-day mass, and ~ 60% upward adjustments based on the luminous stellar mass. The 
"extra" metals are explained by a larger fraction of stars in the supernova-progenitor mass 
range, and a larger ratio of mass in stars formed to present-day stellar mass. 

Successful top-heavy models are characterized by a large "prompt" fraction of SNIa and 
prodigious galactic winds, and have low lock-up fraction and modest asymmetry between 



8 Effects of relaxing the latter are discussed shortly. 
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stellar and ICM abundance patterns. That is, the extra ICM metals are mostly associated 
with the rapid early star formation epoch where both formation of SNcc and SNIa progeni- 
tors, and delivery of the metals from star forming sites to extragalactic hot gas, are efficiently 
realized. 

Successful bottom-light models are characterized by the standard prompt SNIa fraction 
and more modest (though still substantial) galactic winds, with a larger fraction of ICM 
enrichment occurring during the passive, post-star-formation phase. In these models, the 
implied fraction of mass originally in galaxies is 0.42 and the fraction of the ICM that is 
primordial, f pr im = 0.65 - compared to 0.25 and 0.83 in the canonical model. For the top- 
heavy models these take on more extreme values - each on the order 0.5 - although, due to 
the effects of galactic winds, the star formation efficiencies are not appreciably different. 

The different abundance patterns predicted in bottom-light and top-heavy models is 
reflected in values of Rjcm that are ~twice as high for the former, more co nsistent with - 



though still lower than - the value recently inferred by (IBulbul et al.ll2012l ) for the central 
region in Abell 3112. 



4.4. Other Variations 

4-4- 1- Supernova Yields 

In most cases, varying the SNIa yield set primarily affects the predicted (stellar and 
ICM) Ni/Fe ratios - which are not well-determined in clusters at this time. Exceptions 
ar e yield sets with pa rticularly low (< 0.4 M ) Fe yields, i.e. the C-DEF and C-DDT 



m 



Maeda et al.l (120101 ) . These models cannot self-consistently produce the observed stellar 
abundances and Z Fe j C M = 0.5. Since there is a narrow range of Fe yields in SNcc calcula- 
tions, our results are insensitive to the choice of Z cc - though, in principle, abundance ratios 
among a-elements could carry SNcc yield diagnostic information. Decreasing the SNcc upper 
mass limit, m up , with respect to the IMF upper limit, rrihi, lowers the predicted metallicities 
- though the effect is small unless the upper mass IMF slope is very flat. 

Overall our models are conservative in the sense of maximizing Fe yields by adopting 
high values for the SNIa Fe yield (0.74 M ), and for m hi (150 M ). 



Maoz fe Mannucci 



44.2. SNIa Efficiency 
( 120121 ) estimated that a 1.7 x higher Type la supernova rate per unit 
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Table 3. Abundance Ratios in Models with Zp e ,icM = 0.5 









O/Fe 


Mg/Fe 


Si/Fe 


Ni/Fe 


stars 






2.0 


1.4 


1.6 


2.2 


salpeter 






1.0 


0.68 


0.97 


2.8 


canonical 






1.1 


0.77 


1.0 


2.7 


sl-0 


a - 


= 1.83 


1.9 


1.3 


1.5 


2.2 


sl-2 


aio 


= 1.06 


1.2 


0.85 


1.1 


2.7 


sl-3 


ahi 


= 1.97 


1.7 


1.2 


1.4 


2.4 


(e Ia = 0.13) 


&hi 


= 2.22 


0.85 


0.57 


0.87 


2.9 


(star/gas=0.25) 


c*hi 


= 2.36 


1.0 


0.69 


0.97 


2.8 


sl-4 


a-hi 


= 1.85 


1.9 


1.3 


1.5 


2.3 


m-1 


mio 


= 0.20 


1.2 


0.85 


1.1 


2.7 


m-2 


mio 


= 0.42 


1.2 


0.85 


1.1 


2.7 


m-3 




= 1.25 


1.2 


0.85 


1.1 


2.7 



Note. — Canonical and Salpeter model ratios, and stellar ratios, included 
for comparison purposes; e Ia = 0.13 and star/gas=0.25 variations of model 
sl-3 also included. 



Table 4. Characteristics of Models with Z Fe j C M = 0.5 









r* 


frem 


fla 
Jp 


&GW 


£gal 




v Ia 


rf c 


pec 


a SN 


n ICM 


fprim 


salpeter 






0.27 


0.11 


0.5 


0.5 


0.25 


0.17 


0.0014 


0.0065 


0.50 


2.0 


0.32 


0.86 


canonical 






0.41 


0.17 


0.5 


0.5 


0.25 


0.17 


0.0022 


0.011 


0.39 


2.0 


0.27 


0.83 


sl-0 


a - 


= 1.83 


0.56 


0.24 


0.96 


1.7 


0.60 


0.22 


0.0020 


0.019 


0.16 


1.0 


0.11 


0.44 


sl-2 


&lo 


= 1.06 


0.54 


0.22 


0.5 


0.97 


0.42 


0.21 


0.0028 


0.014 


0.24 


2.0 


0.23 


0.65 


sl-3 


&hi 


= 1.97 


0.54 


0.23 


0.78 


1.5 


0.54 


0.22 


0.0023 


0.018 


0.18 


1.3 


0.13 


0.50 


(e Ia = 0.13) 


Chi 


= 2.22 


0.45 


0.18 


0.33 


0.75 


0.31 


0.17 


0.0039 


0.013 


0.32 


3.1 


0.40 


0.76 


(star/gas=0.25) 


&hi 


= 2.36 


0.38 


0.16 


0.46 


0.33 


0.43 


0.32 


0.0021 


0.0096 


0.46 


2.2 


0.32 


0.72 


sl-4 


&hi 


= 1.85 


0.55 


0.24 


0.93 


1.7 


0.59 


0.22 


0.0021 


0.019 


0.17 


1.1 


0.11 


0.45 


m-1 


mio 


= 0.20 


0.54 


0.22 


0.5 


0.97 


0.42 


0.21 


0.0028 


0.014 


0.24 


2.0 


0.23 


0.65 


m-2 


mio 


= 0.43 


0.54 


0.22 


0.5 


0.98 


0.42 


0.21 


0.0028 


0.014 


0.23 


2.0 


0.23 


0.64 


m-3 




= 1.25 


0.53 


0.22 


0.5 


0.98 


0.42 


0.21 


0.0029 


0.014 


0.23 


2.0 


0.23 


0.64 



Note. — Canonical and Salpeter models included for comparison purposes; e Ia = 0.13 and star/gas=0.25 variations of model sl-3 also 
included. 
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mass of star formed is implied by the ICM Fe abundances for otherwise standard assumptions. 
We construct a scenario along these lines by considering an increase in e Ia - the fraction 
of 3 — 8 M Q that explode as SNIa - from 0.076 to 0.13 {f*/fwM an d stellar abundances 
constrained to match the standard values; Section 4.1). Indeed we find Z Fe j C M = 0.5 
for an IMF that otherwise (i.e., in terms of IMF, star formation efficiency, SNcc lock-up 
fraction and rate, and ICM primordial fraction) approximates the standard model (Tables 
3 and 4). Naturally, the increase in efficiency of formation of SNIa progenitors results in a 
large asymmetry between the stellar and ICM abundance patterns (a SN = 3.1), as reflected 
in the more nearly solar \a/Fe icm ratio and consistent with ICM abundance patterns 
(Rfg M = 0.4; iBulbul et all (J201J)) - see Tables 3 and 4, and Figures 9 and 10 that display 
results for such a variation of model si- 3. If this is the correct explanation for the observed 
ICM abundances, one must seek an astrophysical explanation for boosting e Ia in rapidly 
star-forming systems. 



44.3. Star-to-ICM Ratio 

As briefly discussed in Section 2.2, the present-day star-to-ICM ratio may exceed our 
standard value of 11%, e.g. due to an unaccounted-for ICL fraction or underestimate of the 
stellar mass-to-light ratio. If we increase f*/ fwM from 0.11 to 0.25, a generally satisfactory 
resolution of the cluster elemental abundance paradox is achieved^ - a level of Fe enrichment 
and abundance pattern {Rjcm = 0.32) consistent with observations is attained for an IMF 
and other parameters in line with expected values - with the notable exception of the increase 
in star efficiency to e s f = 32% - see Tables 3 and 4. The results of this variation on model 
si- 3 are also displayed Figures 9 and 10. It is worth pointing out at this juncture that we 
defined e s f as the fraction of cluster baryons that form stars; in our models the fraction of 
galactic baryons that form stars (where "galaxies" are defined as locations where stars form 
and eject mass into the ICM) is e s f/e ga i = (1 + 5gu0 _1 - 75% for this /*/ fwM = 0.25 model, 
but also 2/3 for the standard model. 



5. Summary and Discussion 

Star formation, with a canonical IMF and standard efficiency in producing SNIa, that 
builds up a stellar population comprising ~ 10% of the current overall cluster baryon content 
falls short by a factor of > 2 of explaining a typical rich cluster half-solar ICM Fe abundance 



9 Stellar abundances are unchanged. 
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(Sections 2.3-2.4). This is the case even if predicted ICM abundances are enhanced by 
increasing the efficiency at which metals are ejected from galaxies (and where, as a result, 
the overall abundance in stars is significantly below solar), unless the conversion efficiency 
of cluster baryons into stars is also increased well above 10% (Section 3). 

Section 4 (and Appendix A) constructs and utilizes a phenomenological model for the 
evolution of an old, simple stellar population to quantify the changes in the IMF shape 
(high and low mass slopes, break mass) from its standard form required to bring the ICM 
metallicity and cluster stars into concordance in the sense that they be consistent with the 
same parent star formation history. The necessary departure may be either in the "bottom- 
light" or "top-heavy" sense, with the former tentatively preferred based on better agreement 
with observed ICM abundance patterns and on a higher primordial ICM fraction. It is 
further demonstrated that if a standard IMF is to be preserved, a boost in the efficiency 
of forming stars from gas well beyond that consistent with a gas-to-star ratio of 10, and/or 
of producing SNIa progenitor systems, is required. These calculations are conservative in 
the sense of maximizing the enrichment of the stellar population through the choice of SN 
parameters (e.g., SNIa Fe yields, the upper mass limit for SNcc). 

Stars born in cluster potential wells (or those of their progenitors) must be responsible 
for the high level of enrichment measured in the ICM; however, there is increasing tension 
between this truism and the parsimonious assumption that the stars in the generally old 
populations studied optically emerged from the same formation sites during the same epochs. 
Quantifying this tension, and bolstering the case against the universality of star formation 
are the two primary implications of this study. In the remainder of this section we elaborate 
on these themes. 



5.1. The ICM-Enriching Stellar Population as Distinct from that Observed in 

Elliptical Galaxies 



In some cases the departure from the canonical IMF is modest - a shift of a few tenths 
in the slope over some mass range or an increase in the mass at which the slope steep- 
ens from 0.5 to 1.25 M . However, optical determinations of the IMF from kin ematic and 



population studies in elliptical galaxies are trending in the opposit 
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kinematic evidence that indicates a larger mass-to-light ratio in massive ellipticals than ex- 
pected based on a standard IMF is consistent with either an IMF that is top-heavy and so 
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produces more stellar remnants, or one that is bottom heavy and produces more unevolved 
low mass stars. However population synthesis modeling of elliptical galaxy spectra favor the 
latter, and a consensus appears to be emerging for an IMF that is bottom-heavy in elliptical 
galaxies with central velocity dispersions > 150 km s" 1 where most of the present-day stellar 
mass in galaxy clusters reside, being at least as steep as a Salpeter IMF if characterized by 
a single slope (and steeper still at the highest galaxy masses). 

The chasm between the amount of metals expected to be produced from a stellar pop- 
ulation with such a steep IMF, and the observed level of cluster enrichment is illustrated in 
the plots of Fe abundance in the ICM versus cluster star-to-gas ratio for three distinct IMFs 
in Figure 11. Standard yield sets and values of the parameters Sow (0.5), fp a (0.5), and e Ia 
(0.076) are assumed (Section 4.1)0 Also plotted are the overall averaged Fe abundances for 
the total cluster baryons; these are independent of the detailed galaxy evolution parameters 
Saw an d fl a - As might easily be inferred from previous considerations, it is clear that the 
model with Salpeter IMF requires an excessively large gas-to-star ratio, and that an IMF 
as steep as ~ m -3 05 unequivocally falls short by more than an order of magnitude of 
producing the required amount of metals. 

The most straightforward explanation for reconciling the steep IMFs based on optical 
spectroscopic studies of elliptical galaxies and the relatively flat IMFs needed to produce the 
cluster metals is to reject the conventional wisdom that the stellar populations in ellipticals 
that dominate the cluster stellar mass are primarily responsible for ICM enrichment. Such a 
decoupling begs the question of the origin and present-day whereabouts of the ICM-enriching 
stars and motivates consideration of scenarios with pre-enrichment (that would presumably 
be accompanied by pre-heating) in protocluster environments by a currently inconspicuous 
stellar population. 

However, an important caveat with respect to the optical spectroscopic studies is their 
general confinement to regions well inside the half-light radius as well as systems at low red- 
shift. Given the emerging paradigm of the multi-stage, inside-out formation/ assembly /growth 
of ellipticals by multiple mechanisms (star-forming major mergers, "dry" minor mergers, 
and cold and hot gas accretion; e.g., Conselice et al. 2012, Patel et al. 2012), spatial 
gradients and temporal evoluti on in properties of ellip tical galaxy stellar populations such 
as the IMF is to be expected ( La Barbera et allboid ). With the current dearth of global 



constraints, as well as degeneracies betw een the inferred dark matter content and the IMF 
( jWegner et al.ll2012t iTortora et al.ll2012al) a nd possible systematic errors resulting from non- 



solar abundance ratios ( jFerreras et al. 



20121 ). an IMF in cluster ellipticals that is flatter than 



10 



That is, the solutions are no longer constrained to match the standard stellar abundances. 
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currently inferred in the core when integrated over space and ti me is plausible ( jPavel 12008 
Worthey. Ingermann. fc Served 1201 It iNarayanan &: Pavel l2012bl ) . 



Recent arguments for a reconsideration of bimodal star formatio n - hints for an IMF 
that is top-heavy in older, but bottom- heavy in younger, star clusters (jZaritsky et al.ll20!2l ). 
or top-heavy in denser environments (IMarks et al.l 120121 ) - support this. The evolution in 
the cluster star formation environment plausibly leads to an IMF in elliptical galaxies (and 
their progenitors) that transition from one initially weighted towards the mass range that 
includes SNcc and (promp t) SNIa progenitors to on e especially conducive to low mass star 
formation at later times (INarayanan &: Pavel 2012alJbr). We note th at the high metallicities 



seen in the gas in high redshift quasar hosts (IPietrich et al.ll2003al Jbl) also indicate an early 



epoch of rapid star formation that efficiently produces SN progenitors and is accompanied 
by powerful galactic outflows (IPi Matteo et al.ll2004t IWang et al.l 120021 ). 



5.2. The Enriching Stellar Population as Distinct from that Observed in our 

Galaxy 

The hypothesis that star formation is universal is refuted by analysis of the level and 
pattern of ICM elemental abundances. The star formation characteristics of the stellar 
population responsible for these metals must depart from that studied locally in one or more 
of the following ways: (1) engender a higher fraction of high mass stars, (2) more efficiently 
form stars from gas, (3) more efficiently produce SNIa progenitor systems. In addition, we 
saw in the previous subsection that the IMF if the enriching population is distinct from that 
recently inferred in the central regions of elliptical galaxy. Arguments for (1) were presented 
above. We now examine the feasibility, and implications, of hypotheses (2) and (3). 

The true star-to-gas ratio (and implied star formation efficiency) remains uncertain, 
with stellar masses difficult to estimate given low surface brightness extended light and 
the likelihood of multiple stellar populations that complicate the conversion from measured 
light in some aperture to total stellar mass (jMunshi et al. 2012 ). as well as the uncertain 



ICL contribution. Both bottom-heavy as now being inferred in elliptical cores, and top- 
heavy/bottom-light as required by ICM enrichment, IMFs may result in upward revisions in 
mass-to-light ratios. A global value of f*/ fwM > 0.25, even for rich clusters, does not seem 
to be excluded by observations at this time. The star formation efficiency corresponding 
to /*/ fioM > 0.25, excluding the primordial portion of the ICM that does not engage 
in star formation (see Section 4.4.3, above), is {5e ga i(l — f*)} -1 , where e ga i (first defined in 
Appendix A) is the fraction of baryons initially in star-forming structures. This quantity has 
an absolute minimum of 0.2, is > 0.3 for any reasonable value of the mass return fraction 
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r* > 1/3, and > 0.5 for r* > 1/3 and e 9a / < 0.6. These considerations indicate that a 
substantial fraction of protocluster gas was in the form of dense star-forming protogalaxies 
or pre-galactic fragments. There are clearly profound implications for such a high stellar 
fraction and star formation efficiency for e valuating the magnitu de - and perhaps even 
the reality - of the "overcooling" problem ( McCarthy et al.ll201ll ) and the physics of star 
formation quenching as it pertains to the cluster environment, as well as for the precision in 
using cluster gas fractions - that must be converted to baryon fractions using a correction 
for stellar content - to constrain the cosmological world model (lAUen et al.ll2008l ). 



Given the uncertainty in the nature and possible diversity of SNIa progenitors, and the 
difficulties in reproducing observed rates (e.g., Toonen et al. 2012, Quimby et al. 2012, and 
references therein), the feasibility of an efficiency of SNIa progenitor formation in galaxy 
clusters that exceeds the standard e Ia = 0.076 is not easily evaluated, but cannot be sum- 
marily dismissed. Recent work in this area provides hints, on the one hand, of a downward 

vmnciAn m -f- Vi r\ rvl /-vV~i • i 1 r\c^ i c\-i- r~i r~\4~ />i . V\i i +■ /-vn \- \~\ r\ e~\\- ~\r\ nr r-\4~ o l *rln /~i v -t rc\ Inn l vi rvriln -\r~-t r /-~t1 nnf nrc 
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2012; 


Maoz et al. 


2012; 



that produces additional stars in the 3 
value of 77 . 



8 M range, and an increase in e Ia , may boost the 



5.3. Future Directions 



Progress in resolving the cluster elemental abundance paradox will proceed, in parallel, 
along theoretical and obse rvational lines as foll ows. Since there is data on the spatial dis- 
tribution of stars and gas ( jBattaglia et al.ll2012l ). and on the evolution of the Fe abundance 
( jBaldi et al.ll2012l ). we are extending our modeling to multi-zone and time-dependent treat- 
ments - with particular attention to possible mechanisms of pre-enrichment and predictions 
for cluster SN (and 7-ray burst rates; see below) as a function of redshift and - that fur- 
ther constrain enrichment scenarios. We will also extend our investigation to galaxy groups, 
including fossil groups. 

SN surveys are attaining better statistics, particularly at high redshift, and are sharpen- 
ing the accuracy of SN rates, delay-time distributions, and environmental dependencies. Op- 
tical spectroscopic studies are improving both observationally, and in terms of the complexity 
of the stellar population models used to interpret them. X-ray studies of elliptical galaxy 
i nterstellar and circumstellar gas provide additional probes of elliptical galaxy evolution 
( ILoewenstein &; Davisll2010l |2012| ). Future improvements in measuring cluster abundance 
patterns beyond Fe, and in abundance and abundance pattern gradients and time-variation 
are crucial. 
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Finally, we note that many of the mechanisms suggested here for explaining the level 
of ICM enrichment - pre-enrichment by massive stars, efficient and rapid conversion of stars 
to gas, an IMF skewed to high masses - wo uld suggest that the protocluster environment 



is a fertile one f or producing 7-ray bursts ( ILloyd-Ronninget al.l |2002| ; IWang &: Pail 12011 



Elliot et al.ll2012l ). a suggestion we are following up on. 



5.4. Concluding Remarks 

The goal of this work is to quantify the requirements for the stellar population responsi- 
ble for injecting metals into the ICM, and evaluate the feasibility that the stars we see today 
originate from the same source. One is driven to conclude that there is a profound diver- 
gence between the ICM-enriching population and that in the ensemble of elliptical galaxies 
based on standard assumptions about, and recent optical spectroscopic population studies 
of, the latter. This is inferred from the number of SN progenitors needed for the former 
and that expected in ellipticals based on their integrated light and apparent bottom heavy 
IMF, implying the existence of a distinct "hidden" stellar source of metals that may or not 
inhabit the same space as these galactic stars at the same time. While the modeling here 
is basic, the conclusion depend mostly on simple accounting of metals and unlikely to be 
altered in more sophisticated treatments. And although the rich galaxy clusters we con- 
sider represent an extreme environment, there are broader implications f or ellipticals, since 



mass i s a much stronger determinant of their formation than environment (iGriitzbauch et al. 



2011al Jbh. However in it is in the ICM where these phenomena are embedded and remain 



accessible, given the dominance by elliptical galaxies of cluster light, and the closed-box 
nature of these deepest of potential wells. 

We present compelling evidence for a diversity of star formation in terms of some combi- 
nation of efficiency, IMF, and ability to produce SNIa progenitors. Implications to be further 
explored include possible impacts on using cluster baryon fractions to constrain cosmology, 
converting stellar light to mass, and treating star formation and pre- heating/feedback - and 
evaluating overcooling - in semi-analytic models of galaxy formation. 

Occam's razor is violated in rich galaxy clusters - although metals are made in stars and 
most of the stars we observed are in elliptical galaxies, this stellar population as currently 
understood is evidently not responsible for producing the metals in the ICM. Moreover, the 
nature of the star formation that did produce these metals is clearly very different from that 
we are most familar with, as well as that recently inferred in elliptical galaxies. 



The author wishes to acknowledge useful (and enjoyable) discussions with Richard 
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Mushotzky and Esra Bulbul. This paper is dedicated to my late father, Jerry Loewenstein, 
for reasons that cannot be articulated. 



A. A Simple Model for the Composite Chemical Evolution of Cluster Galaxies 



We approximate the stellar population responsible for enriching the ICM as originating 
from a single, brief, and early star formation episode. As discussed in Section 4, we adopt 



Kroupa 


(2001 


); 


Kroupa et al. 



(120121 ) extending from mi to and normalized so that J m h% dm<j)(m) = 1. Thus <p — kf(x) 
for x = m/mi , where 



fix) 



X 



n 



■3i 



X% <C X ^ 



1,3; 



(Al) 



x\ = c\ = 1, Cj = Cj-iXj J x 4 = rrihi/mio, and k is determined from the normalization 

condition. Thus, in addition to the lower and upper mass limits, three slopes (oti,Qt2, and 
a^) and two break-masses (m,2 = rrii X2 and m 3 = rai xz) must be specified. 



For the canonical IMF (IKroupa et al.ll2012l ) that we adopt as default, (m/ c , m 2 , m 3 , rrihi) = 
(0.07,0.5,1.0,150), where all masses are in M© and («i, a2,«3) = (1.3,2.3,2.3). We consider 
IMFs where we vary either rrii (over [0.01,0.5]), m 3 ([0.5,8]); or «i, a 2 , or a 3 (all over 
[0.3,oo)). 



The mass return fraction for intermediate mass stars is given by 

f™ cc dm<j)(m) Am(ffl) 



f™'" dm(f){m)m 



(A2) 



where m to = 0.9 M Q is the main sequence turn-off mass, m cc = 8 M Q is the lower mass limit 
for SNcc progenitors that we use to delineate "intermediate" and "high" mass stars, and 



where 



Am(m) = m- m rem , wd (m), 



fTlrem wd = 0.394 + 0. 199m 



is the white dwarf remnant mass (IKalirai et al.l 120081 ) . Similarly, for high mass stars 

hi dm(j)(m)Am(m) 



f™^ 1 dm<f)(m)m 



(A3) 
(A4) 

(A5) 
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where Am(m) is now based on an averaged remnant mass (see below), 

Am(m) = m- m r em,SNcc(m). (A6) 

The specific number of SNcc explosions per star formed is 

f mup dmd)(m) 

cc J m c .c ^ v 1 



J^ hl dm<f)(m)m 



(A7) 



A simple chemical evolution model for the composite stellar population in clusters is 
constructed and adopted. We use this to calculate the total mass return from high mass stars 
as the metallicity of the progenitor stellar population is built up, and to connect the ICM 
enrichment parameters with the astrophysics of the formation of cluster galaxies. The model 
is appropriate for stellar populations where conversion of gas to stars is relatively rapid and 
efficient, and so may be applied to cluster galaxies where ellipticals dominate the stellar 
mass and star formation is accelerated in general due to the high primordial overdensity. 
As such, galaxy eyo lution is divided into two epochs: active and passive, and three phases 



( ILoewensteinl 120061 ) : star- forming gas ("ISM"), stars, and non-star- forming gas ("ICM"). 
Note that any hot halo gas - relatively insignificant in mass compared to the "true" ICM for 
rich clusters - is subsumed under the ICM category. In the active phase, all star formation 
and SNcc explosions occur and all the initial (ISM) mass in galaxies is consumed by star 
formation or ejected by galactic winds. In the passive phase, stellar mass return and delayed 
SNIa continue to enrich the ICM. 

The cluster as a whole is treated as a closed box, with mass and metal exchange among 
the phases and metal production by t he stellar componen t . We model the active phase 



essentially following the prescription of iQian fc Wasserburgi (120121 ) for the case of no infall. 



Mass return is neglected in their approach, and we correct for this in the passive phase. 
Respectively, the evolution equations for the mass in stars, ISM, and ICM are as follows: 

dM *=M SF , (A8) 

M SF - M GW} (A9) 
= M GW . (A10) 



dt 

dMisM _ 
dt 

dMjcM 



dt 

It is assumed that the rate of outflow is proportional to the rate of star formation that is, in 
turn, proportional to the ISM mass: Msf = ^sfM IS m, Mgw = SgwM S f = ^sfSgwM IS m, 
so that the solutions to equations (A8)-(A10) are 



M ISM = M GALO e~ xt , (All) 
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dfiJSM 




dt 




M SF( 


dt 


M* 


dfiJCM 


_ M GW 


dt 


Micm 



M* = Mqalq ( 1 + (Wr 1 ^ - e- A '), (A12) 

and 

M /CM = M ICM0 + M GAL0 5 GW (l + 5 GW )-\l - e- xt ), (A13) 

where A = \sf(1 + Sqw), and initial conditions correspond to masses of M G alo in star 
forming gas (presumably, mostly in galaxies) and Micmo in the ICM, and M* = 0. 

The evolution equations of the corresponding "ith" element mass fractions are as follows: 

#f ■ < A14 > 

Mism 

i,ISM ~ (A15) 

and 

-(fiJSM ~ fi,ICM)- (A16) 

Mi t SN (the only source term in the set of equations) is the nucleosynthetic production of the 
ith element 

,sn — 8i t sNMsF, (A17) 

6 hSN = rTfe/f > + f p a V Ia vl\ (A18) 

where, as previously defined (see Section 2), yj a and (yf c ) are the yields per SNIa and SNcc, 
rj cc and r] Ia the numbers of SNcc and SNIa explosions per star formed; and, f^ a is the SNIa 
fraction considered "prompt" in the sense that they occur during the star formation epoch 
(not necessarily part of a distinct prompt SNIa mode). 

Analytic solutions for the metal mass fractions are as follows: 

fi,ISM = 8i,SNKt (A19) 

h, = <W(1 + Sow)' 1 { , t \ (A20) 

1 — e At 

and 

kicM = <W<W(1 + Saw)- 2 A f ° AL ° [1 - e- xt (l + At)] , (A21) 

assuming negligible pre-enrichment of any phase. 

Masses and metallicities at the end of the active phase are assigned according to the 
t — 7- oo limit of equations (A10)-(A12),(A19)-(A21), following the presumption that most 



ii 



Only elements primarily synthesized by massive stars and SNIa are considered here. 



-30 - 



star formation in clusters occurs on a timescale much shorter than the current cluster age, 
and then adjusted for the ensuing passive injection of stellar mass return and "delayed" 
SNIa. The final stellar mass and abundances (including remnants) are, therefore, given by 

M* = M GAL0 (1 - r„)(l + 8 GW )-\ = (1 + 5 GW )- 1 (r ] cc {yD + / p W a ), (A22) 

and the final ICM mass and abundances (elemental mass fractions) by 

M ICM = M ICM0 + M GAL0 (1 + 5 GW )-\5 GW + n) (A23) 

and 

fijCM = /o(l + ^)- 1 [l + ^ + /o(^ + r*)]- 1 {(W + r*)r ? cc ( 2/ n 

+ [(1 + Saw) ~ (1 - r,)#W } (A24) 

where the total mass return fraction is r* = r% 1 + r* m . For the overall baryon metallicity 

f l Mr = e sf (v cc (y'i C }+V Ia y! a ), (A25) 

where the star formation efficiency defined in equation (1) is related to e ga i, the "galaxy 
formation efficiency" (the fraction of baryons initially in star-forming structures) according 
to 

e, f = e ga i{l + 5 GW )-\ (A26) 

where e ga i = / /(l + / ) and f = M gaw /M IC mo- 

This formalism enables us to interpret the supernova lock-up parameters introduced in 
Section 2.4 in the context of the chemical evolution of clusters galaxies and place them on a 
firmer physical footing. Naturally the supernova asymmetry parameter a SN = (fp a ) i.e., 
it is the inverse of the fraction of SNIa that explode during the epoch when cluster stars 
form. Our Sections 3 and 4 default a SN = 2 is consistent with estimates of the prompt 
SNIa fraction (Maoz et al. 2011, and references therein; Grauer & Maoz 2012). The lock-up 
fraction may be expressed as /3 CC = (1 — r*)(l + Sqw)' 1 , and naturally depends both on how 
efficiently stars lose mass, and how efficiently star formation induces galactic winds. 



A.l. Mass Return for Massive Stars 

The active phase chemical evolution model is utilized to calculate the total mass return 
from high mass stars - which can be substantial for top-heavy IMFs. This is motivated by 
the profound impact of metallicity on mass loss in massive stars and, hence, fallback and final 
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remn ant mass (1 Woosley &: Hegerl 120021 ; iNomoto et al.l 120061 ; IZhang et al.l l2008t iFrver et al. 



2012 ). The distribution of forming stars as a function of time and mass (IQian &: Wasserburg 
2012n is 

d 2 N M S f 



(f)(m) 



(A27) 



dmdt ^ v J^ ht dm0(m)m' 

from which it follows from the expression for the metallicity of star forming gas, equation 
(A18), that 



d 2 N 



d 2 N fdf Fe 



dmdf Fe dmdt \ dt J f Fe0 ^ ^ 

where fF e o = <5Fe,sv/(l + $gw)- This distribution is used to calculate the mass return from 
massive stars through numerically computing the average remnant mass: 



Jo d fFe ^d^WT^rem^Nccim, f Fe ) 



m 



rem,SNcc 



io d fpe cr dm 



d 2 N 
dmdfp, 



(A29) 



where the remnant mas s as a function of mass and metallicity, m rem! sNcc( m , fFe), is adapted 
from (jFryer et al.ll2012l ) (delayed explosion scenario) using Fe as a proxy for metallicity. 



A. 2. IMF-dependence of SNIa Rate 

The number of SNIa explosions per star is expected to vary with IMF as follows: 

J m i u a P dm(j)(m) 



f™ up dm^{m)rn 



(A30) 



imio 

where m^" = 8 M Q , m[" = 3 M Q , and e Ia — 0.076 yields the observ ationally estimated 
fraction of 3-8 M stars that explode as SNIa (IMaoz fc Mannuccill2012l ). 
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Fig. 6. — Clockwise from upper left (a)-(d): ICM Fe abundance and Mg/Fe, Ni/Fe, and Si/Fe abundance 
ratios for models with departures, with respect to the standard IMF, in either high- or low-mass slope: 
Aa — cthi — 2.3 for top-heavy/light (sl-1, sl-2) and ai — 1.3 for bottom-light/heavy (sl-3, si- 4 , sl~4), models 
- see text and Table 2 for details. Curves for models sl-3 and sl~4 closely trace each other in the ratio plots. 
The solid and broken horizontal lines, respectively, show the stellar, and standard model ICM, values (Table 
I)- 
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Fig. 7. — Same as figure 6 for departures, with respect to the standard IMF, in lower mass limit or break 
mass - see text and Table 2 for details. 
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Fig. 8. — ICM abundance ratios, with respect to Fe, for Mg - left panel (a), Si - middle panel (b), and 
Ni - right panel (c) - vs. Fe abundance for top-heavy models sl-0, sl-3, sl-4, and sl-5; and, for bottom-light 
models si- 2 and m-3 (all bottom-light models considered essentially follow the same curve). 
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Fig. 9. — ICM Fe abundance and Mg/Fe, Ni/Fe, and Si/Fe abundance ratios versus slope above 0.5 M 
(assuming slope 1.3 for m < 1 M Q ) - i.e., model sl-3, for boosted SNIa progenitor formation efficiency (e Ia ) 
or star-to-ICM ratio (f*/ficAi)- As for all models, stellar abundances (solid horizontal lines) are fixed at 
their standard values (Table 1). Results for a standard IMF (broken horizontal lines), and for model sl-3 
with standard e Ia and f*/fic'M, are reproduced. 
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Fig. 10. — ICM abundance ratios with respect to Fe for Mg - left panel (a), Si - middle panel (b), and Ni 
- right panel (c) - vs. Fe abundance for model sl-3 and its boosted e Ia and star-to-ICM ratio counterparts. 
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Fig. 11. — Fe abundance in the ICM (blue curves) and, overall , in cluster baryons (red curves) for the 
followin g three IMFs: (solid curves) standard (jKroupa et al.ll2012l) . (dotted curves) single slope with (j> ~ 
m~ 2 ' 35 ( Salpeter 19551) . (dashed curves) single slope with <f> ^ 
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